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An overview is given of chiral symmetry restoration at finite temperature and bary- 
ochemical potential. Within hadronic models of the vector correlator its implications for 
low-mass dilepton spectra in ultrarelativistic heavy-ion collisions are discussed. 
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§ ■ 1. INTRODUCTION 

ON 

In the limit of vanishing quark masses QCD exhibits an exact chiral symmetry, in which 
left- and right handed quarks decouple giving rise to conserved vector- and axialvector 
currents. In nature this limit is most relevant for the up- and down quark flavors, where 
masses are small compared to typical hadronic scales. To a somewhat lesser extent it 
also applies to the three-flavor case (Nf = 3) which includes the strange quark. Chiral 
symmetry is fundamental to our understanding of light hadron masses where confinement 
^ ■ seems to play a much lesser role. In the physical vacuum chiral symmetry is spontaneously 
broken. For the light meson spectrum representing elementary quark- ant iquark modes 
(Fig. Ill), this manifests itself in two ways: 

- the appearance of (nearly) massless Goldstone particles (pion, kaon, eta). 

- the absence of parity dublets, i.e. the splitting of scalar- and pseudoscalar, as well 
as vector- and axial vector mesons. 

QCD sum rules relate meson masses to the quark condensate (qq)o, an order parameter 
of spontaneous chiral symmetry breaking. For the p- and ai-meson, for instance, one 
obtains in the large A^-limit [H 



2 



448 3 /- V2 

— 7T a s {qq) 



1/3 

m„ 



2816 , 

7r a. 

27 



™\2 



1/3 

(1) 



1.1. Chiral Symmetry Restoration 

It is expected that chiral symmetry is restored at finite temperature, T and quark 
potential, \i q . The only reliable way to obtain the order of the transition, the critical 
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Figure 1. Low-lying spectrum of selected meson states. 



temperature and the critical chemical potential is ab initio lattice simulations. While, at 
H q = 0, results are available with T c = 150 ± 20 MeV H there are difficulties for finite 
fiq. It has long been known that quenched QCD breaks down at non-zero \i q in that the 
transition takes place at \x q ~ m vr /2 which is unphysical. Within random matrix models 
this failure is now understood as an inherent problem of the quenched approximation || . It 
is not clear at the moment whether an early onset of the transition truly holds in full QCD 
or if it reflects computational difficulties in the method employed for the simulation [[|. 

In spite of these difficulties with lattice QCD, model-independent results can be ob- 
tained at low temperatures and small baryochemical potential, /i, from phenomenology. 
Considering the QCD-free energy 



-T 7 = -T In Zqcd 



(2) 



where Zqcd is the grand-canonical QCD-partition function, the condensate ratio can be 
expressed via the difference 5T = JF — JF as || 
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Here m q denotes the quark masses, which act like the external field in a magnet. For a 
dilute gas of hadrons the ratio can be expressed in the two-flavor case as f| 
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where oc dmh/dm q denotes the Sigma-Commutator of hadron h and p s h its scalar den- 
sity. For low temperatures and small baryon density degenerate nucleons and thermally 
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excited pions give the dominant contribution such that 

*S „ ~ 1 " 4^ (^PliT) + E N p s M) . (5) 



ii- 



*o,t^o (qq) mlfl 



Since = = 69 MeV and = 45 ± 8 MeV are known experimentally, the dilute 

gas limit is well determined. The underlying physical picture is very simple. When- 
ever a hadron is created in the vacuum the (qq) condensate is changed locally since the 
condensate inside a hadron is different from that in the vacuum. 

At n = 0, one can go a step further and derive a rigorous low-temperature expansion 
by means of chiral perturbation theory |7]||. Taking the chiral limit — > one obtains 

(QQh . 1 2 16 3 T f T V ra\ 

- — :— = 1 — x x x In— + ... ; x = —= — (6) 

(qq)o 6 9 a \VsfJ 

where the temperature scale is set by V8f w ~ 250 MeV. While the T 2 -and T 4 -terms are 
model independent, model dependence enters at order T 6 through the regularization scale 
A ~ 470 MeV. The range of validity is restricted to T < 120 MeV mostly because at this 
point heavier mesons start to enter. 

For finite fi the dilute gas expression (|5|) predicts a decrease of the condensate ratio 
which is linear in the vector density (for heavy particles the scalar density becomes equal 
to the vector density). At nuclear saturation, p = 0.16//m 3 , this yields a ~ 35% drop 
and naive extrapolation would indicate chiral restoration at ~ 3po- This clearly cannot 
be trusted, since the equation of state of nuclear matter greatly differs from that of a free 
Fermi gas at such high densities pi. 

1.2. How to Detect Symmetry Restoration? 

The principal tool for observing QCD phase transitions are ultrarelativistic heavy-ion 
collisions (URHIC's). In the study of chiral restoration the quark condensate (qq)^T is 
not directly measurable, however, and one has to look for other observables. It follows 
from chiral symmetry alone that, at the phase boundary, the scalar and pseudo-scalar 
correlators as well as the vector- and axialvector correlators must become identical. In 
principle, the observation of parity mixing can thus serve as a unique signal for chiral 
restoration. The study of phase transitions via correlators is known in condensed matter 
physics, as 'soft-mode' spectroscopy [Q. Here structural phase transitions in crystals are 
analyzed by measuring the frequency of low-lying lattice phonons ('soft modes'). A good 
example is SrTiO^ 0. Above a critical a temperature T c = 106 K this crystal has cubical 
symmetry with a three-fold degenerate 'i?25-mode'. Upon cooling this mode becomes 
soft with vanishing frequency at T c . At this point the crystal undergoes a structural 
transition to a phase of tetragonal symmetry. As a consequence the three-fold degenerate 
soft mode splits into two modes one of which is two-fold degenerate. Upon further cooling 
the frequencies of these modes as well as their splitting increase. Obviously the low-lying 
phonon spectrum contains the information of the underlying symmetry change during a 
phase transition. 

The direct QCD- analog is the behavior of the scalar- and pseudoscalar correlators. 
As manifest in the linear a-model they are the carriers of chiral symmetry. As T c is 
approached from above, the scalar (a)- and pseudoscalar (n) modes become soft and split 
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below T c with the pion remaining massless in the chiral limit. The measurement in heavy- 
ion collisions of the scalar- and pseudoscalar correlators is difficult, however, since pions 
suffer strong final-state interactions. A better probe is the vector correlator since it couples 
to photons or dileptons which only undergo electromagnetic final-state interaction. Since 
the vector correlator in the vacuum is saturated by narrow 'resonances', the p, u and 4>- 
mesons, the objective is then to study the spectral changes of vector mesons as a function 
of p and T. There are, in principle, two possibilities. As the scalar and pseudoscalar 
modes, also the vector modes could become 'soft' at T c , giving rise to 'dropping masses'. 
This is the hypothesis of 'Brown-Rho scaling' [plj and would be a natural consequence of a 



direct relationship between the masses and the chiral condensate, as found in the vacuum 
(see eq. (|1])). Brown-Rho scaling also explains qualitatively the rapid increase in entropy 
density across the phase boundary, seen in lattice QCD The second possibility is that 
the vector mesons remain massive at T c , becoming degenerate with their axial partners. 
Support for this possibility comes from the fact that the masses of the p- and ai-meson 
change as T 4 [0 which is difficult to reconcile from a direct link between masses and the 
chiral condensate in the medium. Recall that the condensate ratio, in the chiral limit, 
decreases as T 2 , for small T. 

2. HADRONIC DESCRIPTIONS OF THE VECTOR CORRELATOR 

2.1. General properties 

The vector correlator in the vacuum is defined as the current-current correlation func- 
tion 

KM = * J ^^(01^(^(0)10), (7) 

where T denotes the time-ordered product and j M is the electromagnetic current. For 

three flavors the current can be decomposed in terms of the quantum numbers of the 
relevant vector mesons as 

^ = 3^+3^+3^- (8) 
The corresponding quark content is 

Jl p) = ~{u^u-d^d), (9) 

jf ] = \{u llJt u + d ltl d), (fO) 

J'f = -\{si,s). (ff) 
Current conservation implies a transverse tensor structure 



IWg) = ( g,u - qj f\ I% 2 ), I% 2 ) = ^IV(g). (12) 



The imaginary part of II(g 2 ) is proportional to the cross section for e + e — > hadrons 
. a(e + e~ —> hadrons) 127r . , . . 

= -^r^ = Imn s > 13 

cr(e+e — > p + p ) s 
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where \fs is the total cm. energy of the lepton pair. The vector mesons appear as reso- 
nances in the different hadronic channels carrying the respective flavor quantum numbers 
and the vector dominance model (VDM) quantitatively describes the resonant contribu- 
tion to R(s) [|T3J. 

In a hadronic medium the vector correlator looses its simple vacuum Lorentz structure 
due to a preferred frame. In the grand canonical ensemble one has 

n^ T) (9o,g) = i J d 4 xe^TT^-^ T j,(x)UO))/Z (14) 

with separate dependence on energy and three-momentum. In addition, the longitudi- 
nal and transverse components of II^^ now become distinct. The same holds for the 
propagators 

D^P = ^-{mlf-Ii^P)-\ (15) 

2.2. Theoretical approaches 

At SpS energies the phase space in the final state is dominated by mesons (mostly 



pions) with a meson/baryon ratio of 5-7 |L4|. In calculating modifications of vector 



meson properties through hadronic interactions, an obvious starting point is therefore a 
pure pion gas. For the p meson which in the VDM couples dominantly to two-pion states 
this implies a modification of the pion loop through the heat bath. The resulting thermal 



broadening has been calculated e.g. for on-shell p mesons |15j or within the hidden gauge 
approach [|l^] and found to be rather small. Accounting for off-shell propagation in e.g. 
resonant 7rp-scattering through the ai(1260) meson gives rise to a 7rai-loop in the p-meson 
self energy. Especially for energies above 1 GeV this contribution might be substantial, 
depending, however, crucially on the choice of effective Lagrangian []17|]. 



In spite of the scarcity of baryons in the hot hadron gas, they have a significant impact 
on the spectral properties of vector mesons, largely because of strong meson-baryon cou- 
pling. Several approaches have been put forward to determine these effects reliably. One 
approach starts from a low-density expansion of the vector correlator ]T3| , |r8|| . To leading 
order in density the vector correlator takes the form 

n^' T) (*) = U% + (x) + (ir\Tjp{x)j v (p)W)P*(F) + (N\T Jfl (x)UO)\N)p N (fi) + ■■ (16) 

and involves the pion- and nucleon Compton amplitudes. In the work of Steele et al. Jl8| 
these are obtained from data in a 'chiral reduction formalism' and automatically contain 
mixing of the vacuum vector- and axialvector correlators in the pion sector, as implied by 
chiral symmetry. The nucleon Compton amplitude is constrained from 7iV photoabsorp- 
tion data and the nucleon polarizabilities. In ref. |13| the T = limit of the low-density 
expansion (|16|) is considered. Then only the nucleon contributes and its Compton tensor 
is evaluated in the VDM combined with chiral SU(3) dynamics, based on an effective 
meson-baryon Lagrangian. In the heavy-baryon formalism, the vector meson-baryon am- 
plitude is calculated for p, u>- and meson scattering and significant modifications of the 
spectral properties are found (Fig. ||). Especially the p meson suffers substantial broad- 
ening, as density increases. The position of the 'pole mass' is not affected, however. The 
in-medium uj meson also broadens and exhibits a downward shift of the 'pole mass'. The 
(f> meson spectral function is mostly affected in the sub-threshold region. 
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Figure 2. Modifications of the vector-meson spectrum [13]. Panels (a), (b) and (c) depict 
the p-, to- and (p meson, respectively, The short-dashed lines denote the vacuum values, 
while the dashed (full) lines indicate the results at p /2 (p )- 



A second approach, which has been applied so far only for the p meson, starts from 
the well-known observation that pion propagation in the nucleus is strongly modified. A 
wealth of elastic 7r-nucleus scattering data has provided detailed understanding of the rel- 



evant physical mechanisms fL9 |. The dominant contributions originate from non-resonant 
and resonant ttN scattering. Thus the pion self energy acquires large contributions from 
iV-hole and A-hole loops, giving rise to a momentum-softening of the pion dispersion rela- 
tion. In the VDM it is therefore natural to account for this effect by replacing the vacuum 
two-pion loop with in-medium pions |20T-f22]j . Gauge invariance is ensured by including 



appropriate vertex corrections. To lowest order in nucleon density this approach repre- 
sents a pion cloud model for the nucleon Compton amplitude which coincides with that 



of ref. and ref. jnj (in the latter case there are additional box diagrams). In addi- 
tion to leading-order contributions in p^ there naturally emerge higher orders in density, 
most notably p 2 N terms. These correspond to two-nucleon processes of meson-exchange 
character and iVW- and iVA-bremsstrahlung contributions. Besides the medium modi- 
fications caused by dressing the intermediate 2-pion states, direct interactions of the p 
meson with surrounding nucleons in the gas have to be considered. While elastic scatter- 
ing, pN —>■ N —>■ pN, is kinematically strongly disfavored such a kinematic suppression 
will be much less pronounced with increasing energy of the resonance in the intermedi- 
ate state. Indeed, there are at least two well-established resonances in the particle data 
table |23| which strongly couple to the pN decay channel, namely the iV(1720) and the 
A(1905). This led to the suggestion p3| to consider p-like particle-hole excitations of the 



Table 1 

Baryon resonances with appreciable decay width into pN. 



B* 


I(JP) 


I> [MeV] 


iV(1520) 


1/2(3/2") 


24 


A(1620) 


3/2(1/2-) 


22.5 


A (1700) 


3/2(3/2-) 


45 


iV(1720) 


1/2(3/2+) 


105 


A(1905) 


3/2(3/2+) 


210 



7 



type piV(1720)iv"~ 1 and pA(1905)iV _1 . In a more complete description other resonances 
with appreciable pN widths have been included [f25|-|27jj. These are listed in Tab. 1. The 
most simple version of the VDM overestimates the B* — > N'-f branching fractions when 
using the hadronic coupling constants deduced from the B* — * Np partial widths. This is 
corrected for by employing an improved version of the VDM [p8| , which allows to adjust 
the B*N , ~f coupling independently. 

Combining the effects of pionic modifications and resonant pN scattering, the resulting 
p- meson spectral functions IhlD p = 1/3 [D P L + 2Dj) are displayed in Fig. ^| 




Figure 3. The p-meson spectral function [25,27] in the vacuum (left panel), at pn = 
p , T = (middle panel) and at Pn + Pa = l-4po ; T = 166 MeV (right panel); notice that 
the vertical scale in the middle and right panel is reduced by a factor of 2. 

One observes significant broadening especially at small q with a low-mass shoulder 
which originates from intermediate AiV -1 7r states and resonant iV(1520)iV~ 1 excitations. 
As q increases the shoulder moves towards the M = line, which can be understood from 
simple kinematics. As seen from the right panel, temperature has little influence on the 
energy-momentum dependence and only provides an additional overall broadening. 

2.3. Constraints from photoabsorption 

An obvious constraint for models of p-meson propagation is photoabsorption in nuclei, 
for which a wealth of data exist over a wide range of energies. Real photons correspond 
to the M = line in the middle part of Fig. |3|. Within the model discussed above, the 
total photoabsorption cross section per nucleon can be calculated straightforwardly. Tak- 
ing the low-density limit, pn — > 0, only terms linear in density contribute, representing 
the absorption process on a single nucleon. Two contributions have to be distinguished: 
(i) coupling of the photon to the virtual pion cloud of the nucleon, incorporated in the 
VDM via pion dressing; this contribution describes the background contribution the ab- 
sorption process, as mentioned above; (ii) resonant contributions via the intermediate B* 
resonances listed in Tab. 1. 

Adjusting the model parameters to optimally reproduce the data |27J yields results 
displayed in the left panel of Fig. |j. Photoabsorption on nuclei can be reproduced with 
similar quality. It is noteworthy that, in the nucleus, strength below m v is obtained, 
which originates from two-nucleon processes via meson-exchange currents and is nothing 
but the 'quasi-deuteron tail' of the giant dipole resonance. 
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Figure 4. The photo absorption spectrum on the proton (left panel) and on nuclei (right 
panel) as obtained in ref. [27]. The dashed lines indicate the non-resonant background 
contributions. The data were taken from refs. [29-33] . 



2.4. Comparison with Dilepton data 

Constraining hadronic models via photoabsorption [IS, 27] gives confidence in extrapo- 
lations to the time-like region of dilepton production. For the p meson the dilepton rate 
is obtained as 



dN„ 



+i+i- 



« 2 K) 4 f p (qo;T) 



d 4 xd 4 q 



3n 3 g 2 



M 2 



(17) 



In the of model |25j most of the important processes discussed above have been included 
such that the p-meson propagator becomes 



f)P npTTT p)pNB* | Tjpnai , y)pKKx 



;is) 



and contains aside from the baryonic contributions also the most important parts from 
the pion/kaon gas. 

To compare the theoretical rates with data, the space-time history of the heavy-ion 
collision has to be specified. There are several possibilities for modeling the collision. 
The most naive approach is a simple 'fireball' model p3,Z5|, in which initial conditions 
in temperature and hadron abundances are obtained from transport model calculations. 
Assuming local thermal equilibrium as well as chemical equilibrium, the space-time history 
is then determined by a simple cooling curve, T(t), from some initial time, up to the 
'freeze-out time', tf. 

For SpS energies such cooling curves are available from transport theory |35| and can 
be easily parameterized. The time evolution of the hadron abundances is determined by 
chemical equilibrium and agrees well with transport model results. The observed spectrum 
is then obtained by integrating the 'local rate' (|i~7"|) in time, accounting for the detector 
acceptance in addition. Results from [25, 18 1 are displayed in Fig. ||. While the full results 
in the right panel [p5] give a reasonable account of the observed spectra, especially in the 
region below the p, to mass, the results of ref. |18| fail to describe this region. The reason 
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Figure 5. CERES Dilepton Spectra [36] employing the expanding fireball model of ref. [25] 
with in-medium dilepton rates from ref. [18] (full lines in left panels) and ref. [25] (full 
lines in right panels). 



for this discrepancy between the two calculations is currently not completely understood. 

The fireball model is rather crude and does not incorporate detailed flow dynamics. 
A more satisfactory description is provided by hydrodynamical simulations. In this case 
the 'local rates' (|TT|) refer to a given fluid cell and can be directly implemented. So far, 
however, no results are available. 

A third way is to supply transport calculations with rates obtained from in-medium 
vector-meson propagation. For the p meson this has been implemented in simulations us- 
ing the HSD model for the transport f3"7|| . Results have been obtained for both the 
CERES e + e~-measurement in Pb+Au collisions (left panel of Fig. as well as the 
HELIOS-3 ^^"-measurement of S+W collisions (right panel of Fig. §). in both cases, the 
in-medium rates give an improvement in the comparison between theory and experiment 
and the overall description is quite satisfactory. 
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Figure 6. SpS Dilepton Spectra [36,38] employing the HSD transport model [37] using 
free (upper panel) and in-medium [25,27] (lower panel) p-meson spectral functions. 

3. SUMMARY 

Recent theoretical efforts in understanding the nature of chiral symmetry restoration 
at finite temperature and baryochemical potential and their implications for low-mass 
dilepton production in URHIC's have been discussed. While the onset of chiral sym- 
metry restoration is well understood in the dilute gas limit, the approach towards the 
phase boundary within hadronic models remains a theoretical challenge. Here symmetry- 
conserving nonperturbative methods are called for. A unique signal for chiral symmetry 
restoration in URHIC's would be the observation of degenerate vector- and axialvector 
correlators. Such a mixing strictly follows from chiral symmetry but is difficult to detect. 
Only the vector correlator is accessible via electromagnetic probes. 

The application of hadronic models to assess in-medium properties of vector mesons 
and their consequences for dilepton spectra have reached a degree of maturity such that 
realistic comparisons with data become meaningful. While some approaches inherently 
involve constraints from chiral symmetry [ 13|Jl8 , others lack an obvious connection to 



chiral symmetry emphasizing, however, input from hadronic phenomenology [p4 , 25 ,P7 . 
When supplemented by models for the space-time history of the collision, quantitative 
agreement between theory and experiment seems to emerge fl2"5 , j3~7| . Further calculations 
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addressing more exclusive observables remain to be confronted with experimental data on 
a quantitative level. 
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